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Abstract. We revisit the index leverage effect, that can be decomposed into a volatility effect and a 
correlation effect. We investigate the latter using a matrix regression analysis, that we call 'Principal 
Regression Analysis' (PRA) and for which we provide some analytical (using Random Matrix Theory) 
and numerical benchmarks. We find that downward index trends increase the average correlation 
between stocks (as measured by the most negative eigenvalue of the conditional correlation matrix), 
and makes the market mode more uniform. Upward trends, on the other hand, also increase the 
average correlation between stocks but rotates the corresponding market mode away from uniformity. 
There are two time scales associated to these effects, a short one on the order of a month (20 trading 
days), and a longer time scale on the order of a year. We also find indications of a leverage effect for 
sectorial correlations as well, which reveals itself in the second and third mode of the PRA. 



1. Introduction 

Among the best known stylized facts of financial markets lies the so-called "leverage effect" [131 [H 
[131 HE1 E] ! a name coined by Black to describe the negative correlation between past price returns 
and future realized volatilities in stock markets [5] . Q It is indeed well documented that negative price 
returns induce increased future volatilities, an effect responsible for the observed skew on the implied 
volatility smile in stock option markets (see e.g. 3, 4, 11 ). 

However, the association, made by Black, with a true leverage effect (i.e. that when the value of 
a stock goes down its debt to equity ratio increases, thereby making the company riskier and more 
volatile), is probably misleading. In particular, the amplitude of the leverage correlation for indices is 
noticeably stronger than for individual stocks, which even sounds paradoxical when the index return is 
by definition the average of individual stock returns! The volatility of an index in fact reflects both the 
volatility of underlying single stocks and the average correlation between these stocks. The increased 
leverage effect for indices must therefore mean that both these quantities are sensitive to a downward 
move of the market. 

The aim of the present paper is to investigate more specifically this "correlation leverage effect" , 
and make precise the common lore according to which correlations "jump to one" in crisis periods (see 
[12j EU HH EQ] for early studies of the time evolution of the correlations in financial markets). Similar 
studies have appeared recently. In [T] , a careful study of the average correlation between stock returns 
during contemporaneous upward/downward trends of the market index has confirmed that correlations 
are indeed stronger when the market goes down [7] . Our analyses confirm and make more precise these 
results, first by extending them to different markets, and second by devising and exploiting a new tool 
to investigate conditional correlations, that we call "principal regression analysis" (PRA). The idea 
here is to regress the instantaneous correlation matrix on the value of the index return (or any other 
conditioning variable) . While the intercept of the regression gives the average correlation matrix, the 
regression slopes define a second symmetric (but not definite positive) matrix that can be diagonalized, 
leading to modes (eigenvectors) of sensitivity to the conditioning variable(s). The interpretation of 
these eigenvectors is particularly transparent when they coincide with those of the correlation matrix 
itself. The corresponding eigenvalues quantify how the whole correlation structure of stock returns is 
affected by the conditioning variable. The nice point about the PRA is that Random Matrix Theory 
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^While this effect holds for most markets in developed economies, Tenenbaum et al. [22] report that the situation 
appears to be different for markets in developing countries. 
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(RMT) provides, as for standard PCA, a useful guide to decide whether or not these sensitivity modes 
are statistically meaningful (for a review on RMT, see [TU]). When the conditioning variable is the 
past values of the index return, the conclusion of PRA is that the dominant mode is the market 
mode, associated to a negative eigenvalue, indeed corresponding to a correlation leverage effect. We 
characterize the temporal decay of this effect. Upon separating positive and negative index returns, we 
furthermore find that the correlation leverage effect is strongly asymmetric: whereas negative returns 
increase both the volatility of the underlying stocks and the average correlation between stocks, positive 
returns have weaker influence on these quantities (see Fig. 6 below). We furthermore find indications 
of a leverage effect for sectorial correlations as well, which reveals itself in the second and third modes 
of the PRA. 



2. Data, notations and definitions 

We have considered 6 pools of stocks corresponding to 6 major stock indices: SP500, BE500, Nikkei, 
FTSE, CAC 40 and DAX. We analyze the daily returns in a time period spanning from 01/01/2000 
to 04/26/2010. Stocks are labelled by a — 1, . . . , N (where N depends on the market), and days by 
t = 1,...,T (where T = 2594). Time average will be denoted by (.). The return of stock a between 
the close of day t — 1 and the close of day t is denoted as r/ a (i) . We in fact understand ij a (t) as the 
demeaned return over the whole time period T. We define an inverse volatility weighted index return 
at time t as: 

'(*) = 4$>«(t), = W 

a— 1 

where a a is the average volatility of the stock a over the whole time period: 

«£:=-£*.(*)'. (2) 

4=1 

We will further define the average instantaneous stock volatility a(t) at time t as: 

1 N 

*(*)*•■= nY,**® 2 (3) 

a=l 

while the average instantaneous correlation between all pairs of stocks p(t) is defined as: 

p ( t )-= 1 v MOM! (4) 

The average over time of the above two quantities will be denoted as cr 2 , and po- 

The squared index return I(t) 2 is a rough proxy for the instantaneous index volatility. Using the 
above definitions and the fact that N is large, it is easy to check that: 

i(t) 2 «p(fMt) 2 + 0(i), (5) 

showing that both the average stock volatility and the average correlation contribute to the index 
volatility. It is therefore natural to decompose the full index leverage effect in two contributions: one 
coming from the dependence of the average stock volatility on the past returns of the index, and a 
second one describing the average correlation. We thus define a full leverage correlation function Cj{t): 

Ci{t) - — (jw) — ' ( } 

and two partial leverage correlation functions: 

_ (I(t-r)a(tr) (I(t-r)p(t)) 
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All the above leverage correlation functions are normalized to be the regression slope of the corre- 
sponding observables on the past value of the index return, for example: 

p{t) = p + C p (r)I{t-T) + e(t,T), (8) 

where e(t,r) is some noise. (Remember that by construction, I(t) has zero mean.) 
In the limit of weak correlations, the two effects are additive and one should find: 

C t {t)^ p Q L a {r) + (xlC p {T), (9) 

eliciting the contribution of the average stock volatility and of the average correlation to the full leverage 
correlation. The second term is responsible for the enhanced leverage effect for indices compared to 
single stocks. 




FIGURE 1. Left: normalized leverage correlation functions po£ a (T), <7q£ p (t), and an 
exponential fits with two scales (dotted lines). Right: Full leverage function £/(t) 
and comparison with an additive model (dotted line). 



3. Index leverage effect: A simple empirical analysis 

As a first stab at understanding the index leverage effect, we plot in Fig. [T] the normalized partial 
leverage correlation functions, po£ a (T), <Tq£ p (t), together with the full leverage £i(t). In these plots, 
the data is averaged over the four indices, SP500, BE500, Nikkei and FTSE. From this figure, we draw 
the following conclusions: 

• (a) the two contributions to the index leverage are of the same order of magnitude. In partic- 
ular, the correlation leverage is significant and confirms the conclusions of Refs. PQCj. 

• (b) the correlation effect is stronger at short times but decays faster than the volatility effect; 
a two time scale exponential fit of these two contributions in the range r € [1,250] (in days) 
indeed leads to 

cjIC p {t) ks -0.053 exp(-r/18) - 0.005 exp(-r/350); (10) 
PoC a {r) w -0.02 exp(-r/14) -0.02 exp(-r/280), (11) 
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• (c) a test of Eq. Q with the sum of the above two fitted exponentials reproduces satisfac- 
torily the full leverage effect, although the latter is underestimated at short times, when the 
correlations cease to be small enough for Eq. ^ to be accurate. 

In fact, one can test directly whether linear regressions such as Eq. Q above make sense or not, by 
averaging all values of pit) corresponding to a given value of lit — 1) within some range. The resulting 
graphs are shown in Fig. [2J both for p and for a 2 . One sees that whereas a linear regression for p 
makes sense for I(t — 1) < 0, there is in fact perhaps a small positive slope for I(t — 1) > 0. For 
o" 2 , the graph looks even more symmetric, reflecting the presence of volatility correlations on top of 
(assymetric) leverage correlations. 
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FIGURE 2. Dependence of the average correlation p(t) and the average single stock 
volatility cr 2 (i) on the index return the previous day, I(t — 1). The result is obtained 
as an average over all 6 indices: SP500, BE500, Nikkei, FTSE, CAC40 and DAX, but 
the qualitative effects are robust and appear on each markets individually. These plots 
suggest that a quadratic I 2 (t — 1) term should be included to the linear regressions. 
The printed error bars are the average of the error bars obtained for each of the 6 
indices. 



4. A MORE PRECISE TOOL: THE "PRINCIPAL REGRESSION ANALYSIS" 

The above analysis, although interesting, is oversimplified, because the structure of inter-stock 
correlations is described by a full correlation matrix C and not by a single number p, that only 
captures the average correlations. In order to characterize the way the correlation matrix depends on 
the past value of the index (or on any other conditioning variable) , we propose the following: consider 
a given pair of stocks, a, j3, and regress the product of normalized returns rj a (t)rjp(t) on the past value 
of the index return, i.e. write: 

rja{t)rj{s{t) :=C at p + D a ,p(T)I(jt-T)+e atfi (t,T). (12) 

Since I(t) has zero mean, the intercept of the regression is exactly the empirical Pearson estimate of 
the correlation matrix. The regression slopes D a ^(r) define another N x N symmetric matrix D(t), 
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which encodes the full information about the dependence of the correlations on past returns. More 
precisely, the regression leads to the following empirical determination of D(t): 

1 T 

(I 2 )D aiP (r) = — — -r £ %(t)f} P (t)I(t-T). (13) 

^ T > t=T+l 

The aim of this section is first to discuss the information contained in D(r), in particular its eigenvalues 
and eigenvectors, and second to use results from Random Matrix Theory to assess how meaningful this 
information is when the length of the sample, T, is not very large compared to the number of stocks 
N. Finally, we describe our empirical results on D(r), in particular its most negative eigenvalue and 
eigenvectors. 

4.1. Interpretation. Define C(J) to be the correlation matrix conditioned to a certain past value of 
I, by: 

C(J) = C + 7D. (14) 

The interpretation of the matrix D is particularly simple when it commutes with the correlation 
matrix C, i.e. when the eigenvectors of D are the same as those of C. In this case, the eigenvectors 
of C(J) are exactly the same as those of C, whereas the eigenvalues Afc(I) are shifted asj^] 

A fc (/)=A fc (0)+/<« fc |D|i; fc ), (15) 

where Afc(O) are the eigenvalues of C and \vk) are the associated eigenvectors (in quantum mechanics 
notations). When D does not commute with C, the structure of the eigenvectors themselves is impacted 
by the conditioning variable. If DJ is small enough, standard first order perturbation theory gives 



back Eq. ( 15 ) for the eigenvalues and: 



A ' Afc — AP 

for the eigenvectors of the matrix C{I). 

As we will find below, the eigenvector corresponding to the most negative eigenvalue of D turns out 
to be very close to the first eigenvector of C (i.e. the so-called market mode, whereas all other 

eigenvalues are significantly smaller. In this case, the top eigenvalue of C is to a good approximation 
given by: 

Ai(7)« Ai+Z/ii, (17) 

where p\ is the most negative eigenvalue of D. Since Ai can be used to define the average correlation 
between stocks through Ai := Np, the meaning of p± is similar to, but more precise than, the correlation 
leverage function C p defined above. 

More generally, when D and C do not commute, one expects the "correlation leverage" to rotate 
the top eigenvector away from the market mode \v\) . The common lore is indeed that when markets go 
down, all stocks "move together" , meaning that the top eigenvector should rotate towards the uniform 
vector |e) = (1/%/iV, 1/y/N, . . . , 1/VN). The cosine of the angle between \v\) and |e) is given by the 



scalar product (e|t>i), that one can compute using perturbation theory. Eq. (16) above. Assuming 



further that the top eigenvalue of C is much larger than all the others (Ai 3> A^i), one finds: 

(e\ Vl (I)) « (e\ Vl ) + f [<e|DM - (vx]^} (e\ Vl }] . (18) 

A measure of how strongly the top eigenvector moves towards |e) is therefore provided by the quantity 
A, defined as: 

A=^[(e\B\v 1 )-(v 1 \B\v 1 )(e\v 1 )}. (19) 

A negative A means that the instantaneous market mode is closer to the uniform mode |e) when the 
index goes down, since (e|ui(/)) — (e\vi) = I A > 0. 

2 Note that the dependence on the lag r is implied in the following formulas. 



6 



PIERRE-ALAIN REIGNERON, ROMAIN ALLEZ, AND JEAN-PHILIPPE BOUCHAUD 



4.2. Results from Random Matrix Theory. When N is large, the simultaneous determination 



- using Eq. (13) above - of the N(N + l)/2 different elements of D from the NT data points is 
problematic, exactly in the same way the correlation matrix C is hard to measure. We thus need to 
provide a benchmark to compare the empirical results obtained with the noise level of the benchmark 
case. This will enable to separate significant effect from noise level arising from the dimensionality 
problem. Let £ be a random variable which will play the role of the conditioning variable (the past 
values of index returns in our context) and let x a ,a = 1,...,N be a gaussian vector of covariance 
matrix C which should be seen as instantaneous stock returns. The x a will be supposed to have 
mean and unit variance, so that C is the correlation matrix of the gaussian vector (x±, . . . ,Xn)- 

We begin by the case C = I. Suppose, in addition, that there is no correlations whatsoever between 
the conditioning variable £ and the correlation x a xp, and that one forms a matrix D from: 

1 T 

(^A^y £>»(*)*/»(*)£(*)• (20) 
t=l 

In the limit T — > oo for finite N one should find that all the elements of the matrix D are zero, 
and therefore all its eigenvalues are zero as well. For finite T, however, the matrix D will have a 
set of non trivial eigenvalues. Random Matrix Theory offers a way to compute the statistics of these 
eigenvalues when TV and T are both large, with a fixed ratio q = N/T. The result depends both on the 
eigenvalue spectrum of the matrix C and, perhaps surprisingly, on the probability distribution of the 
conditionning variable, P(£). The simplest, albeit unrealistic case for applications in finance, is when 
C is the identity matrix, i.e. there is no correlations between the rj. In this case, using the theory 
of Free Random Matrices [23], one finds that the empirical eigenvalue spectrum of D, px(fi), is the 
solution of the following set of equations, in the limit where e goes to zero: [TO] E] 

m = r^^ 2 + ! dgP(0 h f r ~!f, R ] p „ (2D 

Gr + k Pi J I 1 ~ kGr) + KQ^Pi) 

qe 



Pi ( - m , ^ + / CT m_^, 2, ~* ) . (22) 



where Gr is the real part of the resolvent. One can check that in the Unfit q — > 0, and using the fact 
that £ has zero mean, the above equations boil down to: 

- — = (i-ie-t p 1 (ii)=6(p), (23) 

Gr - inpi 

i.e. all eigenvalues are zero, as they indeed should when T 3> N. 

The case of an arbitrary correlation matrix C can also be solved completely using the above result 
on pi and the so-called ^-transform of the eigenvalue spectrum [23], noting that the eigenvalues of 
D are the same as those of the product C x Di, where Di is a random matrix with eigenvalue 
spectrum p\ (ji) . The resulting equation can in principle be solved numerically for any value of q and 
for an arbitrary correlation matrix C. The resulting theoretical eigenvalue spectrum for the matrix D, 
assuming no correlation between the conditioning variable £ and the instantaneous correlation x a X/3, 



can be compared to the empirical spectrum obtained from data using Eq. ( 13 ). Any difference between 
the two spectra can be interpreted as resulting from a true correlation with the conditioning variable. 
In the null-hypothesis case, it is also clear that the quantity A defined by: 



e|D|«i)-(«i|D|t;i)<e|wi) . (24) 



must be zero when averaged over £, x a . One can compute its variance, which is found to be: 



<e|c V; <eM V ). (25) 
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For large T, the central limit theorem ensures that A becomes Gaussian with the above variance. This 
result will be used below to assess whether the empirical value of A (defined above) is meaningful or 
not. 




-40 -30 -20 -10 



Figure 3. Main figure: empirical spectrum of D for the BE500 index (in red), com- 
pared to the null- hypothesis case (in blue). For the latter case, we have generated 
1000 random samples, ranked the eigenvalues and averaged each of them separately. 
The leftmost blue peak therefore corresponds to the average value of the most negative 
eigenvalue. Insets: cumulative distributions of the most negative eigenvalue /ii and of 
the scalar product S — (wi\vi). 



4.3. Numerical simulations. In practice, however, we found it more convenient to use direct nu- 
merical simulations rather than the above exact results. In principle, these results below could be 
obtained using the mathematical formalism above, but the effort required to solve numerically the 
equations above is larger than the one needed to make direct simulations. We measure the null- 
hypothesis spectrum of D by choosing to be a Gaussian random variable of zero mean and unit 
variance, completely independent of the true returns r) a (t), which we then diagonalize. The cumulative 
distribution of the largest negative eigenvalue in the null- hypothesis is shown in the inset. The average 
position of the most negative eigenvalue of D in the null-hypothesis case is found to be Jli s» —4.8. 
The average position of the second and third most negative eigenvalues in the null-hypothesis case will 
be denoted by fa and fa. 

We have also measured the distribution of the scalar product S = (wi\vi) between the corresponding 
top eigenvector \wi) and the top eigenvector of C, \v%). We find that even in the case where is an 
independent random variable, the top eigenvector of D is in fact strongly correlated with \v\), with an 
average scalar product equal to S = 0.68 for the correlation matrix of the returns of the BE500 index. 
We find numerically that P(S < 0.5) w 0.11 and P(S < 0.65) w 0.38 for the BE500 index - see Fig. [3} 
Results for the SP500 are very similar. 



4.4. Comparison with empirical data. In order to reduce the measurement noise and compare 
with the above numerical simulations, we have estimated D(r) using Eq. (13) with "Gaussianized" 
empirical index returns, obtained by first ranking the true index return from most negative to most 
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positive, defining the rank of day t, k(t). The Gaussianized index return Ic(t) is then obtained as 
$ -1 (fc(i)/T), where $ is the error function. 




i i i i i i ^ i i i i f 

50 100 150 200 250 50 100 150 200 250 



FIGURE 4. Left: Largest negative eigenvalue /ii(r) of the lagged regression ma- 
trix D(t). The double exponential fit (dotted line) is given by : /ii(r) = fif — 
26.6 exp(— r/f 1) — f 7.1 cxp(— t/200), where we fix the value of using the numer- 
ical results of the previous section: fxf = fix ~ —4.8, since we expect that for large 
t, all correlations are lost. Right: Evolution of the scalar product S(t) = (v%\wi(t)) 
as a function of r. The horizontal dashed line corresponds to the mean of the scalar 
product S in the null-hypothesis case. The data corresponds to the BE500 index, but 
the results for the SP500 are very similar. 



We show in Fig. [4] the evolution of //i(r), the largest (in absolute value) eigenvalue of D(t) as a 
function of r. We find that fi\ is negative, corresponding to the correlation leverage effect (see Eq. (17)). 
Comparing with the null-hypothesis case, we find that (J,i(t) remains significant at the 1% confidence 
level up to t ~ 240. When fitting /ii(t) with an exponential function with two scales that saturates at 
the noise level ~\x\ determined above, we find \i\(r) =fx\ — 26.6 exp(— r/11) — 17.1 exp(— r/200). This 
reveals two time scales; a rather short one close to the one determined directly from C p (t) above (see 
Fig. [lj, and a much longer time scale on the order of a year, showing that the effect of market drops on 
the correlation is long lasting. The scalar product S(r) — (wi(t)\vi) between the top eigenvectors of 
D(t) and C globally exceeds 0.8 in the whole range r € [1,240], whereas the null- hypothesis average 
value is S = 0.68. 

We have also studied the second (/^(t)) and third (^(r)) eigenvalues of D(r) as a function of r, 
which are both negative and clearly beyond the noise level, and are found to decay with very similar 
time scales: a month and a year (see Fig. |Eb. The corresponding eigenvectors are found to be mostly 
within the subspace spanned by the second and third eigenvectors of C. The financial interpretation of 
these eigenvalues is of an increased sectorial correlation when the market drops on top of an increase of 
the market correlations. Therefore, all idiosyncratic effects disappear upon market drops, while global 
factors become dominant. 



4.5. Separating negative &; positive returns. As Fig.[2]explicitely shows, the correlation depends 
on past index returns in a non-linear way. In fact, both negative and positive returns increase the 
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FIGURE 5. Left: Second eigenvalue ^(t) of the lagged regression matrix D(r). 
The exponential fit (dotted line) is given by: ^(t) — jj-2 — 1.3 exp(— r/14.4) — 
2.3 exp(— t/364). Right: Third eigenvalue /x 3 (r) of the lagged regression matrix 
D(t). The exponential fit (dotted line) is given by : /^(r) = JI3 — 1.3 exp(— r/20) — 
1.5exp(— t/420). Direct numerical simulations of the random case lead to JI2 ~ —1-52 
and JI3 « —1.24. The data corresponds to the BE500 index, but the results for the 
SP500 are again very similar. 



correlations, although the effect is stronger for negative returns, which in turn leads to a non-zero 
linear term in the regression of r} a (t)Vi3(t) on I{t ~ T )- A wa Y to capture the parabolic shape seen in 
Fig. [2] would be to extend the above model to: 

V«{t)Vp(t) := C aS + D aJj (T)I(t - r) + E a>0 {r) [I 2 (t - t) - (I 2 )] + e a ^(t), (26) 

defining a new matrix E that captures the symmetric effect of index returns on the correlation matrix. 
An alternative choice, that we adopt below, is to regress separately on negative returns and on positive 
returns: 

rfa(t)rtp(t) := C aJi + D+ p (r) [l+(t - r) - (/+)] 5 {I{t _ T>0} (27) 
+ D *A T ) [ /_ (* - r ) - < 7 ")] *{iM<0} + £ °At)> (28) 

where / + = max(J, 0), / _ = min(J, 0) and 5 is the Dirac function. With this definition, one can rewrite 
the correlation matrix conditioned to a certain past value of I more precisely separating the effect of 
positive returns and negative returns, as follows: 

C(J) = C + D- [I- - {I-)] 6 {I<0} + D+ [/+ - (/+)] <5 {/>0} . (29) 

Again, in order to reduce the measurement noise, we used "Gaussianized" empirical index returns Jg(^) 
instead of I(t). We apply to D ± (r) the same analysis as above. As anticipated, the top eigenvalue 

of D~ is strongly negative, whereas the top eigenvalue fi^ of D + is positive, but with /if < 
— see Fig.[(3] The projections of \w^) and \w^[) onto \vi) are both very close to unity for small r and 
gradually decay to the noise level as r increases. To check the significancy of our effect, as before, we 
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define a null-hypothesis case, introducing the matrix: 

1 T 

t—i 

where the conditioning variables <f> is independent of the x a (which are standard gaussian variables 
whose correlation matrix is C as above) and distributed as min(£, 0) — (min(£, 0)) where £ is as before 
a standard gaussian variable. We define further the matrix D + exactly as D~ except for the fact that 
the conditioning variable is now distributed as max(£,0) — (max(£, 0)). As above, /Ij - , , will be 
the average positions of the first, second and third most negative eigenvalues of D and ji^ , ju J , jUg" 
will be the average positions of the first, second and third most positive eigenvalues of D + . Those 
values are all computed using numerical simulations. 




1 1 1 1 1 1 — ^ 1 1 1 1 r~ 

50 100 150 200 250 50 100 150 200 250 



FIGURE 6. Left: /^f (r) of the lagged regression matrix D~(r). The exponential fit 
(dotted line) is given by : Hi(t) = Jl^ — 73exp(— r/19) — 41 exp(— r/300). Right: 
Hi(t) of the lagged regression matrix D + (r). The exponential fit (dotted line) is now 
given by : Hi(t) = Ji\ + 10.6 exp(— r/49) + 44exp(— r/200). Note again the presence 
of a long relaxation time on the order of a year. We have used direct numerical 
simulations to obtain Jl^ = Jlf s» —8.3. The data is for the returns of the BE500. 
Again, SP500 yields very similar results. 

We have also studied the rotation parameter for both matrices D ± (r) defined as: 

A ± = i[<e|D ± |« 1 )-(« 1 |D ± |i; 1 )<eh)]. (31) 
Ai 

The results are shown in Fig. [7] In agreement with the common lore, A - is negative, indicating that 
strongly negative index returns (below lead to a more uniform instantaneous market mode. On 

the other hand, A + is found to be negative as well, meaning that while strongly positive returns also 
tend to increase the average correlation between stocks, the instantaneous market mode rotates away 
from the uniform vector |e). The effects w e ar e reporting are statistically significant since the root- 
mean square error on A* (defined as in Eq.(24|) in the null-hypothesis case is found to be ~ 8 • 10 -4 , 



a 



factor 3 to 4 smaller than the amplitude of the empirical values of A*. 
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Figure 7. Plot of the rotation parameters A~ (left) and A + (right) as a function of 
r. The horizontal dashed lines correspond to the root-mean square error on A ± in the 
null-hypothesis case. The data is for BE500; the amplitude of A - and A + are found 
to be roughly a factor 2 larger for the SP500. 



5. Summary & Conclusion 

The aim of this paper was to revisit the index leverage effect, that can be decomposed into a volatility 
effect and a correlation effect. We investigated the latter in great detail using a matrix regression 
analysis, that we called 'Principal Regression Analysis' (PRA) and for which we have provided, using 
Random Matrix Theory and simulations, some analytical and numerical benchmarks. 

Using this refined analysis, we confirm that downward index trends increase the average correlation 
between stocks (as measured by the top eigenvalue of the conditional correlation matrix), which in 
turn explains why the index leverage effect is stronger than for single stocks. Compared to the null- 
hypothesis benchmark, this leverage correlation effect is highly significant (see Fig. [4] and Fig. [6]). We 
also find that large downward trends implies a more uniform future market mode (see Fig. [7j left). 

Upward trends, on the other hand, also increase the average correlation between stocks (see Fig. [6j 
right) but large upward trends rotate the future market mode away from uniformity (see Fig. [7j right). 
All these effects are characterized by two 'memory' time scales: a 'short' one on the order of a month 
and a longer one on the order of a year. The latter long time scale could be related to the fact that the 
market had long cycles of booms and busts within the studied time series, during which the average 
correlation went down and up again. 

We have also studied the correlation leverage effect on intraday data, and we find (results not 
shown) that while the top eigenvalue of the 15 minutes correlation matrix is nearly insensitive to the 
sign of the previous 15 minutes index return, a significant effect emerges when the time scale reaches 
one hour. 

Finally we have found indications of a leverage effect for sectorial correlations as well, which reveals 
itself in the second and third modes of the PRA (see Fig. (51). It would be interesting to analyze 
other conditional correlation matrices using the tools developed in this paper, such as for example 
leader-lagger effects 19, 8, 18J, or the role of other macro variables such as oil, currencies or interest 
rates. 
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